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The statistical properties of energy eigenvalues in the critical regime of the lowest 
Landau band are investigated. The nearest neighbour level spacing distribution 
P(s) and the Mehta quantity Io, from which the scaling exponent v of the cor- 
relation length can be derived, are calculated. While P(s) disagrees with recent 
analytical predictions, v is found to be in agreement with theoretical and experi- 
mental results. 



1 Introduction 

In recent years the investigation of energy level statistics in systems showing 
a localisation-delocalisation transition has made remarkable progress. For the 
metallic side of the transition the spectral fluctuations only depend on the sym- 
metry of the system (orthogonal, unitary or symplectic) and are well described 
by Random Matrix Theory (RMT)ld. On the insulating side, the energy levels 
are uncorrelated and exhibit Poissonian behaviour. At ±hj&, transition point, 
new-ifiritical ensembles were found for the 3d orthogonal atra the 2d symplec- 
tic EKJ, and the 3d unitary case 0. In the quantum Hall system na-j^oiprjlete 
metal-insulator transition is observed, and the situation is unclear ffl't2Hiil. 

This is due to the fact that the localisation length of the electronic wave- 
functions diverges at the centre of the disorder broadened Landau bands. 
Hence, in the thermodynamical limit there is no range_pf extended states but 
only one energy where the states are multifractall§HL3. Therefore, the ques- 
tion is whether there exists a scale-invariant level statistics in finite systems in 
the absence of true extended states. 

As a result of analytical theories EjE_3 the asymptotic behaviour of the 
nearest neighbour level spacing distribution P(s), where the energy difference 
s is measured in units of the mean level spacing A, was proposed to obey 
P(s) oc exp(— cs 2-7 ) with 7 = 1 — -K, where v is the critical exponent and d 
the Euclidean dimension. Fits to the numerical data in 3d orthogonal □, 3d 
unitaryd and 2d symplectic Ua systems did not result in a sufficient agreement. 
Investigations that seem to show a reasonable ntLinij require, however, values 
for v which are different from those commonly accepted. Previously, a simple 
exponential decay of P(s) was proposed a, P(s) cx exp(— ns), where n > 1 is a 
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Figure 1: Energy and size dependence of io at the critical energy in the lowest Landau band 

numerical constant. From the calculations in 3d orthogonal^ and 3d unitary 
systems 13, n ~ 1.9, and for the 2d symplectic, k « 4.0 was obtained In 
our opinion, k should reflect the fractal structure of the critical wavefunctions. 
Therefore, it is very interesting to study the asymptotics of P(s) in the 2d 
QHE situation where no complete localisation-delocalisation transition occurs. 



2 Model and method 

We investigate a 2d Anderson model with the magnetic .field included via a 
phase- factor in the transfer matrix elements V (see Rem for details). The 
diagonal disorder potentials are taken at random with constant probability 
from [-W/2, W/2] and W/V = 2.5. The strength of the magnetic field equals 
1/8 flux quanta per unit cell. We diagonalise the corresponding Hamiltonian 
matrix for finite squares of size L 2 — 32 2 , 64 2 , 128 2 , and 160 2 using a Lanczos 
procedure. To obtain a reliable statistic, up to 3 6 realizations were calculated. 

3 Results and discussion 

To find the energy range of the eigenstates with localisation length larger than 
the system size and to extract the critical exponent v 1 we employ the quantity 
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Figure 2: The critical QHE-P(s) (•) in comparison with the critical 3d unitary (*), the 
Poisson and the GUE result. The inset shows a exp(— ks) decay with k ?s 4.1 

Iq known from RMT0. The calculatiomdoes neither require a fitting procedure 
as the Brody level repulsioivparameterB, nor some cutoff value, as for example 
the integrated A3 statistics Ell. Iq is obtained by integrating over P(s) 

/•OO pOO />OC 

I = / / P(s")ds"ds'ds. (1) 

JO Js Js' 

The value of Iq for the Gaussian unitary ensemble is Iq = 0.590 and for un- 
corrected (localised) levels Iq = 1. Multifractal states are expected to have a 
value larger than 0.59. From our calculations we find Iq = 0.605 ± 0.002 near 
the critical energy independent of the considered system sizes. Fig. ^ displays 
the energy dependence of Iq in the energy range about the critical point. We 
have carefully checked that Iq does not dependent on the width of the energy 
interval from which the eigenvalues are taken. For all system sizes one ob- 
serves a minimum Iq m 0.605 at the critical energy E c /V — —3.33. States 
away from E c become more localised when the system size is enlarged. From 
the scaling behaviour of Iq(E,L) we extract the critical exponent v = 2.4 in 
good agreement with both theoreticalEj and experimentalLfl results. 

The nearest neighbour level spacing distribution P(s) is shown in Fig. ^. It 
clearly deviates from the Wigner surmise although it shows a quadratic small-s 
behaviour as expected for a unitary symmetry. Our result for the lattice model 
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resembles that obtained for a continuum model Due to the large number 
of realizations we are able to extract the large-s behaviour displayed in the 
inset of Fig. ||. As for the othfir critical ensembles, it is not possible to fit the 
tail of -P(s) by the proposalEj P(s) ~ exp(— cs 2-7 ) with v = 2.4, but it can 
be fitted by a simple exponential decay, P(s) cx exp(— ks), with k w 4.1. This 
value is similar to the one found for a 2d system with symplectic symmetryEl. 
In conclusion, there exists a critical level statistics in 2d QHE systems despite 
the absence of a complete metal-insulator transition. 
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